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We present the first study of superconducting nanowires shunted with an external resistor, geared 
towards understanding and controlling coherence and dissipation in nanowires. The dynamics is 
probed by measuring the evolution of the V-I characteristics and the distributions of switching and 
retrapping currents upon varying the shunt resistor and temperature. Theoretical analysis of the 
experiments indicates that as the value of the shunt resistance is decreased, the dynamics turns 
more coherent presumably due to stabilization of phase-slip centers in the wire and furthermore 
the switching current approaches the Bardeen's prediction for equilibrium depairing current. By a 
detailed comparison between theory and experimental, we make headway into identifying regimes 
in which the quasi-one-dimensional wire can effectively be described by a zero-dimensional circuit 
model analogous to the RCSJ (resistively and capacitively shunted Josephson junction) model of 
Stewart and McCumber. Besides its fundamental significance, our study has implications for a range 
of promising technological applications. 

PACS numbers: 74.78.Na, 74.25.Fy, 74.25. Sv, 74.40.+k, 74.50.+r 



I. INTRODUCTION 

The dissipation of the supercurrent in thin supercon- 
ducting wires occurs solely due to Little's phase slips 
[l|. Advances in fabricating ultra-narrow superconduct- 
ing nanowires has greatly boosted the interest in study- 
ing phase slippage in quasi-one-dimensional supercon- 
ductors @. There has been an intense activity to es- 
tablish the existence of quantum phase slips (QPS) re- 
lated to macroscopic quantum tunneling (MQT) [3l-Fllj| 
and to study quantum phase transitions between possi- 
bly superconducting, metallic and insulating phases in 
nanowires [1, [T^ - [l7t. A dissipation-controlled quantum 
phase transition [l8i-25] has been predicted in junctions 
of superconducting nanowires. Recently the importance 
of taking into account Joule- heating caused by dissipative 
phase-slip fluctuations has also been argued and demon- 
strated both theoretically and experimentally [Tol [TTI. |26|- 
[28} . Furthermore, quantum theory shows that the QPS 
rate as well as the quantum phase transition can be 
controlled by an external shunt [23]. Dissipation plays 
an important role in dictating the physics of nanowires. 
Conversely, superconducting nanowires provide an ideal 
prototype for studying the interplay between coherence, 
dissipation and fluctuations. Besides their fundamental 
importance, superconducting nanowires are also ideally 
suited for building superconducting nano-circuitry and as 
devices with potentially important applications, such as 
superconducting qubits and current standards [29|, [30j . 
Thus, even from the technological point of view, it is 
extremely important to fundamentally understand the 
mechanism and role of dissipation in nanowires and to 
find a way of experimentally controlling coherence and 
dissipation. 

It is well established that the environmental dissipa- 
tion of Josephson junction (JJ) can be controlled by 
externally shunting the junction [3l|. This effect has 



been observed in the voltage-current (V-I) characteris- 
tics, which are greatly altered by the amount of dissi- 
pation. The statistics of the switching and retrapping 
behavior in shunted JJs have been investigated in last 
three decades and continues to be actively studied [32h - 
[39j . In general, the retrapping current, which is in- 
versely proportional to the quality factor Q of the circuit, 
is more sensitive to the amount of damping/dissipation 
than the switching current. The shunting is also known 
to control the rate of MQT of the phase variable in 
superconductor-insulator-superconductor (SIS) junctions 
[40l - |47j . The Stewart-McCumber model [H, H9| of resis- 
tively and capacitively shunted Josephson junctions (RC- 
S Js) accurately describes much of the physics of shunted 
JJs [U H(| ■ This model is quite useful since it allows 
the analysis of various fundamental aspects of supercon- 
ducting devices, including chaotic behavior [51j and high- 
frequency microwave responses (52| . The analysis of su- 
perconducting computational circuits also involves use of 
the RCSJ model ||. 

Shunting a superconducting nanowire with a normal 
resistor should have a strong effect on the superconduct- 
ing character of the wire and just as in the case of JJ 
could potentially provide a powerful way to control co- 
herence and dissipation. In spite of its clear importance 
and relevance, the behavior of shunted nanowires has not 
been studied previously both experimentally and theo- 
retically. Here, we present the first study of shunted 
nanowires. It has been inarguably proven for the case of 
unshunted nanowires that going beyond linear response 
is essential to probe the dynamics of (quantum) phase- 
slip fluctuations [rJ EH I2q - l28j . Furthermore, most ap- 
plications would require the wire to be driven out-of- 
cquilibrium, making it doubly important to understand 
how the dynamics of shunted nanowires evolves upon 
shunting. In fact, there is a third equally important mo- 
tivation for such a study. As discussed above, the RCSJ 
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model has been successfully used for JJs and has proven 
to be extremely important. However, a circuit-element 
representation of a superconducting nanowire is currently 
lacking and through this work we want to fill this gap by 
making some concrete advances in that direction. 

The nanowires on which measurements were performed 
in the present work were located in a low-pressure, ther- 
malized helium gas and were fabricated using the molecu- 
lar templating method resulting in suspended nanowires 
[HIH . It had already been demonstrated that these super- 
conducting nanowire show a large hysteresis in the V-I 
characteristics for the unshunted case and that this hys- 
teresis stems from Joule heating and the strong tempera- 
ture dependence of the resistance of the wire [l(J, l26T - l28j | . 
Local Joule heating by phase-slip processes is especially 
important for a long free-standing nanowire because the 
heat generated in the bulk of the wire is not removed 
easily and has to flow away through the ends of the wire. 
Observation of similar physics in a recent study of Alu- 
minium nanowires fabricated using a different method 
fill ] points to the ubiquity and importance of Joule heat- 
ing effects and further underlines how it can be turned 
into an effective probe for quantum phase slips. How- 
ever, the best case scenario will be to be able to have a 
control over the Joule heating. As will be shown in this 
paper heating can indeed be controlled by shunting the 
superconducting nanowire with an external resistance. 

The article is organized as follows. In the next Sec. 
IrEl we briefly describe the sample fabrication and mea- 
surement technique. The experimental results are pre- 
sented in Sec. IIII1 This is followed by our theoretical 
analysis and discussion in Sec. IIVI where we will argue 
that shunting qualitatively changes the behavior of the 
nanowire and present the theoretical results obtained by 
modelling the nanowire. Finally we will end with con- 
cluding remarks in Sec. [V] 



II. SAMPLE FABRICATION AND 
MEASUREMENT TECHNIQUE 

The nanowires presented in this study are fabricated 
using molecular templating 0, 3 ■ Using electron-beam 
lithography and a reactive ion etch, a 100 nm wide trench 
is patterned in the SiN layer of a Si-Si02-SiN substrate. 
The trench is then etched in a 49% solution of hydroflu- 
oric acid to form an undercut to prevent electrical leak- 
age between the electrodes, which are separated by the 
trench 54]. Fluorinated single- walled nanotubes, which 
are insulating, are dissolved in isopropanol and then de- 
posited onto the substrate containing the 100 nm wide 
trench in the SiN layer and then dried with nitrogen gas. 
Randomly, some of the nanotubes cross the trench, cre- 
ating a scaffold for the nanowires to form as the metal 
of choice is deposited on the substrate. The samples are 
then DC sputtered with amorphous Mo 76 Ge24 in a high 
vacuum (~ 10~ 7 Torr base pressure) chamber, thus coat- 
ing the substrate and nanotubes with 12-18 nm of MoGe 



depending on the sample. A scanning electron micro- 
scope (SEM) is then used to image the trench until a 
MoGe coated nanotube (nanowire) is found to be rela- 
tively straight, homogeneous, and coplanar with the elec- 
trodes [2| . An SEM image of one such nanowire is shown 
in the inset of FigJ^fa). Contact pads are formed using 
photo lithography and wet etching in a 3% solution of 
H2O2, which etches MoGe rapidly. 

All of the samples studied in this paper are ~100 nm 
long and are fabricated using MoGe. The thickness of 
each nanowire is controlled by the deposition time in 
the sputtering chamber and by the configuration of nan- 
otubes used as a scaffold. The actual width of each sam- 
ple is measured from the SEM image and found to be 
-15, 12, 10, 15, 8, and 18 nm for samples SI, S2, S3, 
S4, S5, and S6 respectively. Thicker samples show a 
lower normal resistance i? n , higher critical temperature 
T c , higher critical current / c , and slightly higher retrap- 
ping current J r . For example, for samples SI, S2 and 
S3 which have a decreasing thickness, the resistance R n 
in the normal state and critical temperature T c are SI 
(Rn = 1385 ft, T c = 4.607 K), S2 (R n = 1434 ft, T c = 
4.41 K), and S3 (R n = 1696 ft, T c = 3.82 K). All sam- 
ples visually show a similar behavior in the R-T and V-I 
curves. 

The shunt has been added by attaching a commercial 
metal film resistor (ranging from 5 to 200 ft) parallel to 
the sample using silver paste (FigQJb)). The distance 
from the nanowire to the shunt is 1-2 cm for all samples, 
and the shunt resistance is measured as a function of 
temperature down to cryogenic temperatures and found 
to be constant. Measurements are performed in a 4 -He 
or 3 He cryostat equipped with base temperature silver 
paste and copper powder filters and room temperature 
7r-filters. Transport measurements are carried out by 
current biasing the sample through a large resistor (—1 
Mft) and measuring the voltage with a battery-operated 
Stanford SR 560 preamp, using the typical film-inclusive 
four-probe technique as in Fig[TJb) [J]. Resistance vs. 
temperature (R-T) curves are measured by applying a 
small sinusoidal current (~10-100 nA) at a frequency 
of ~12 Hz and measuring the voltage and then doing 
a linear fit to the resulting voltage-current data to ob- 
tain the resistance. The temperature is measured using 
a calibrated Cernox thermometer from LakeShore. V-I 
curves are measured by applying a large sinusoidal cur- 
rent in the range of a few /iA, at a frequency of a few Hz, 
and measuring the voltage simultaneously. The switch- 
ing (retrapping) current has been measured by sweeping 
the current as in the V-I measurement, and recording 
the current at which the voltage-jump (drop) out of the 
superconducting (resistive) state has been the greatest. 



III. EXPERIMENTAL RESULTS 

In this section we present a series of experimental 
results. We begin by presenting the temperature de- 
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pendence of the linear resistance of shunted nanowires 
and use that as a benchmark for characterizing differ- 
ent samples. Next, we go one step further and present 
the measurements of the nonlinear V-I characteristics of 
the current-biased nanowires at temperature of 1.8 K, fo- 
cusing on the evolution of hysteresis upon shunting. In 
the next subsection the experimental data of switching 
and retrapping distributions for different values of shunt 
are presented. Such an analysis provides deeper insights 
into the dynamics of superconducting nanowires as is ev- 
ident from a previous study in the case of unshuntcd 
wires PS IH Although the focus of our work is on 
studying the effect of shunting, in the last subsection we 
present the temperature dependence of the V-I charac- 
teristics and the switching distributions so as to provide 
a comparison with the corresponding measurements for 
the unshunted case that was studied in detail! 10. 27, 281. 



A. Shunt dependence of R(T): characterization of 
the samples 

Study of the linear-response resistance as a function 
of temperature is useful for characterizing the samples 
and establishing a starting point for further investiga- 
tion. Figfjja) shows the temperature dependence of the 
nanowire's resistance using a log-linear scale for various 
values of the shunt resistance i?g. As the temperature is 
lowered below 5.8 K the film becomes superconducting 
while the wire is still resistive because its critical temper- 
ature (T c ) is lower than that of the film. Below T c of the 
wire, as expected there is a measurable resistance due to 
phase slips in the wire. 

We understand the measured resistance vs. tempera- 
ture (R-T) curves using the following arguments. Below 
T c of the nanowire the total sample resistance is a par- 
allel combination of the Rs and the wire resistance, i?w- 
We model the wire resistance with an empirical formula, 
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R W (T) R n R AL (T) 



(3.1) 



where R n is the normal state resistance of the nanowire 
to account for the quasi-particle resistance channel and 
i?AL is the Arrhenius-Little (AL) resistance occurring due 
to thermally activated phase slips (TAPS). The AL re- 
sistance is estimated, following Little's proposal, by as- 
suming that each phase slip creates a normal segment 
on the wire of a size equal to the coherence length and 
for a time interval roughly equal to the inverse attempt 
frequency P, HJ . We note that the Langer- Ambegaokar- 
McCumber-Halperin (LAMH) theory [HI HI of TAPS is 
not valid except very near to T c [251 ] . So we have to use 
the phenomenological AL expression: 

Ral(T) = Rn exp(- ^P), 
where AF(T) = {^ 1 )M(X) 
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FIG. 1: (a) R-T data for sample SI with various shunt val- 
ues. The first sharp drop in resistance at 5.8 K is due to 
the film electrodes going superconducting. The second, grad- 
ual, drop in resistance at lower temperature is the supercon- 
ducting transition of the nanowire, which is much broader 
due to TAPS, (b) Sample schematic. The wire is shown as 
a short vertical line and is shunted by a commercial resistor 
(Rs)- The sample is measured by current biasing the sample 
and extracting the resistance via a four-probe measurement, 
(c) Comparison of R-T curves from (a) for the case of no 
shunt (circles) and a 25 il shunt (squares) with the theoreti- 
cal expression of the total sample resistance given by Eq. (|3.3p 
(solid lines). The known parameters used in each fit are: wire 
length, L=105 nm, normal state resistance, R n — 1385 Q, and 
shunt resistance, Rs = oo and Rs = 25 f2, respectively. The 
fitting parameters used for the fits are T c = 4.607 K for the 
unshunted case and T c = 4.595 K for the 25 SI shunted case, 
and £(0)=7.57 nm for both cases, where £(0) is the dirty limit 
coherence length at zero temperature. 



is the free energy barrier for a phase slip in the zero- 
bias regime [55| . Here H C (T) is the thermodynamic crit- 
ical field, £(T) is the temperature dependant coherence 
length, A is the cross-sectional area of the wire, and fee 
is the Boltzmann constant. The equation for the free 
energy barrier AF(T) can be rewritten to include wire 
parameters more accessible via the experiment as [ItJ 

AnT)=0.83 fcB T c @(^)(l-^) 3/2 , (3.2) 

where L is the length of the nanowire and Rq(— 



4 



h/4e 2 6450 f2) is the quantum resistance. Thus, the 
temperature-dependent total sample resistance is: 

Jhm _£ + £ (lM ^ )r (3 . 3) 

The fits of the total sample resistance by Eq. (|3.3|) are pre- 
sented in FigOTc) f° r the unshunted nanowire (Rg = oo), 
and the case when the nanowire is shunted with 25 f2. 
All the fits are done by using the values of R n , T, L, 
and i?s obtained from the R-T curve and the SEM im- 
age and using T c , and £(0) as fitting parameters. The 
value of the fitting parameters change very slightly as 
the shunt resistance is varied. For instance, in the fit- 
ting presented in Figfljc), the value of T c decreased by 
12 mK for the shunted case compared to the unshunted 
case, which can be accounted for by slight sample change 
during thermal cycling. The agreement between experi- 
mentally measured resistance and the total resistance in 
Eq. p.3[) , as shown in FigQJc), gives evidence that the 
T-dependence of the rate of phase slips does not depend 
on the shunt at relatively high temperatures (T > 4 K in 
this case) and that the observed residual resistance of the 
wire just below T c is due to TAPS in the high tempera- 
ture limit of ~ 4— 5 K. Note also that the R-T curves of 
all the samples presented in this paper are smooth and 
show no extra transitions, and the SEM images confirm 
that the nanowires are homogeneous and well-connected 
to the electrodes. Thus our nanowires are well suited to 
systematically study the effect of shunting. 

B. Shunt dependence of V-I characteristics 

Let us start by considering the unshunted wire (Rs = 
oo). As the bias current I is increased, thermal fluctua- 
tions cause the nanowire to switch from a superconduct- 
ing state into a resistive state before the current reaches 
the critical (equilibrium) depairing current. The current 
at which the wire switches out of the superconducting 
state is called the switching current (I S w)- Once in the 
resistive state, as the current is decreased below some 
critical value of current, the nanowire experiences retrap- 
ping back into the superconducting state. The current 
at which this happens is called the retrapping current 
(7 r ). For unshunted wires, switching is stochastic in na- 
ture, i.e., each new current sweep gives a different value 
for the switching current but the retrapping process is 
non-stochastic. We devote the next subsection for the 
discussion of switching and retrapping distributions and 
their dependence on the shunt resistance. Here we focus 
on the mean values of the switching and retrapping cur- 
rents and show how they evolve upon shunting the wire 
as shown in Figf^l 

FiglUJa) shows the V-I characteristics for different val- 
ues of shunt resistance. As the nanowire is shunted with 
lower values of the shunt resistance, the mean switching 
and retrapping currents are increased while the width of 
the hysteresis is decreased. Additionally, the retrapping 
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FIG. 2: (a) Voltage vs. total current for sample SI at 1.8 K 
corresponding to various values of the shunt resistance, Rs- 
Dashed circles are kinks in the V-I curve occuring for the case 
when the nanowire is shunted with Rs = 25 Q. Inset: SEM 
image of the nanowire for sample SI. (b) Mean switching and 
retrapping currents vs. G s (= 1/Rs) for three samples (SI at 
1.8 K, S2 at 1.5 K, and S3 at 1.6 K). 



current also becomes stochastic (as shown in Fig0]). In 
FigfSJb), the dependence of the mean switching and re- 
trapping currents on the shunt resistance are shown for 
different nanowire samples. The mean switching current 
increases, at a lower rate than the retrapping current, 
and saturates for small values of the shunt resistance 
(FigHJb)) with a decreasing (increasing) shunt resistance 
Rs (conductance Gs)- Similarly, the retrapping current 
increases with decreasing the shunt resistance Rs until R 
finally reaches I sw of the wire. Such behavior is observed 
on all tested samples. As the switching and retrapping 
current coincide, the hysteresis disappears. Nanowires 
with smaller T c start showing this saturation behavior at 
higher shunt values. 

When the nanowire is shunted with a 25 f2 resistor 
or less, kinks in the voltage are also observed (they are 
marked by dashed line circles in Figj2ja) for the case 
where sample SI is shunted by 25 fi. The kinks for the 
10 Q shunting case are not shown, but occur at higher 
current. We relegate the interpretation of these kinks to 
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FIG. 3: Switching distributions vs. bias current for sample SI 
shunted with various _Rs values at 1.8 K. All distributions were 
measured with a sinusoidal current sweep with f — 8 Hz and 
an amplitude of 6.1 /xA. Inset: Shows the standard deviation 
CTsw of switching current distribution as a function of G s (= 
1/ife). 
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C. Shunt dependence of switching and retrapping 
distributions 

As mentioned in the previous subsection, the un- 
shunted nanowire undergoes stochastic switching as the 
bias current is increased. We plot the switching distri- 
butions vs. current in Fig |3] for different values of exter- 
nal shunt resistance to see how it evolves upon shunting. 
Lowering the value of the shunt resistance has the effect 
of narrowing the width, increasing the height, and shift- 
ing the distribution to higher currents. As can be seen 
from the plots, the full width of the distribution at half 
maximum (FWHM) changed from 100 nA to 12 nA due 
to shunting with a 25 57 resistor. The asymmetric shape 
in the distribution for larger shunts changes to a rather 
symmetric shape with lower shunts. 

The retrapping current also shows a dramatic change 
from deterministic values to stochastic values when 
shunted with 75 fi or less. The bottom part of Fig HI is a 
typical retrapping histogram for an unshunted wireflOj. 
Here, the standard deviation of the retrapping current 
is 1.51 nA, which is the noise limit of our experimental 
setup. So, this small distribution of retrapping current is 
just due to the instrumental noise, and can be reduced by 
decreasing the instrumental noise and the spacing in be- 
tween bias-current points. Thus, retrapping in unshunted 
nanowire always occurs at the same current, i.e. the tran- 
sition is deterministic. However, when the wire is exter- 
nally shunted, a retrapping distribution is observed, with 
its width being much larger than the experimental setup 



FIG. 4: Top: Retrapping distributions vs. bias cur- 
rent for sample SI shunted with various Rs values. The 
mean and standard deviation of the retrapping data are: 
3.683, 4.166, 4.925 /jA and 7.304, 6.089, and 5.799 nA, respec- 
tively for Rs = 75, 54.5, 25 Q. Bottom: Typical retrapping 
histogram vs. I for the case where the nanowire is unshunted. 
The standard deviation of the retrapping current of the un- 
shunted wire matches with the experimental current noise. 

noise and independent of the bias-current spacing. 

In the top of FigfJ] the retrapping current distribu- 
tions for sample SI are shown for the nanowire shunted 
with different values of external resistances. The width 
of the distribution is slightly sensitive to the value of 
shunt resistance, but the mean value of the retrapping 
current changes considerably. When shunted with 75 51 
for instance, the standard deviation of the retrapping 
current increases above the experimental setup noise to 
7.3 nA (from 1.5 nA for the unshunted case under the 
same conditions). Interestingly the width of the retrap- 
ping and switching current distributions for 25 57 shunt 
is almost similar. The retrapping distributions for the 
shunted nanowire are asymmetric in contrast to the un- 
shunted case in which the distribution is symmetric as 
can be seen in the bottom of FigfJJ 

D. Temperature dependence: shunted vs. 
unshunted nanowires 

In this subsection we discuss the temperature evolution 
of the dynamics of shunted nanowires. In FigfSfa), the 
mean value of the switching current is plotted at various 
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the equilibrium depairing current, Bardeen's prediction 
for the temperature dependence of the equilibrium 
critical (depairing) current is compared to the tempera- 
ture dependence of the measured switching current for 
nanowires shunted with small resistances as shown in 
FigEl The Bardeen's equation, derived from BCS theory, 
is given by: 



7 c (T) = / c0 (l-(^) 



3/2 



(3.4) 



where I c q is the critical current at zero temperature. Ex- 
cellent agreement is found with the experimental data 
over a wide temperature interval suggesting that the 
shunt has driven the switching current close to the de- 
pairing current. In these fits, the temperature is known 
while the critical current at zero temperature I c q and 
the critical temperature T c are used as fitting parame- 
ters. Close agreement is found between the theoretical 
prediction for the depairing current at zero temperature: 
I c0 = 92 f iA(LT c )/R n £(0) [Hj, which is derived from BCS 
and Ginzburg-Landau theory, and the value of i c o used in 
the Bardeen fit. Here L and £(0) are in nm, T c in K , and 
R n in SI. Using the fitting parameters from the Rt(T) 
fit of Eq. (|3.3l) presented in FigQTc), IcO has a theoretical 
value of 5.48 /iA, while J c o used in fitting to the Bardeen 
formula in Figj5]has a value of 5.88 /iA. Thus, excellent 
agreement is found between the theoretical and experi- 
mental value for 7 c o. A value of 4.2 K for T c is used to 
fit the temperature dependence of the switching current 
with the Bardeen formula, while the Rt(T) fit using the 
AL model predicts the value of T c to be 4.717 K. This 
difference can be accounted for by sample oxidation and 
thermal cycling between measurements. 



FIG. 5: (a) Mean J sw vs. T for samples SI and S4 shunted 
with 5 SI and 10 £7 respectively. The solid lines are fits to the 
Bardeen's prediction fEq |3.4|) for the temperature dependence 
of the depairing current, 7 C . The fitting parameters used are: 
ico = 5.88, 5.53 /iA and T c = 4.2, 3.9 K for samples SI and 
S4, respectively, (b) Distribution of 7 SW vs. T for sample S5 
unshunted (red) and shunted with 30 Q (black). Each curve 
contains approximately 20,000 points. The corresponding fit 
to the Bardeen prediction (green) is presented for the shunted 
wire. Here I c0 = 1.415 /iA and T c = 2.62 K. 
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temperatures for sample SI shunted with 5 SI and sam- 
ple S4 shunted with 10 SI. In FigJSJb), a distribution of 
switching currents is plotted as a function of temperature 
for sample S5 when it is shunted with a 30 SI resistor and 
compared with the unshunted case. As the temperature 
is reduced, the switching current for all samples increases 
and begins to show signs of saturation below 1 K. The 
behavior of the switching current as a function of tem- 
perature in the distribution measurement in FigJSfb) is 
similar to that of FigJSJa) except that in Fig|5jb) the 
fluctuation in the switching current is also displayed. 

To check the proximity of the switching current to 
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FIG. 6: The standard deviation <t sw vs. T for the switching 
events of sample S6 shunted with 10 fl, where for each point 
a was calculated using data sets of 10000 points. 

Finally in Fig 111 we plot the measured temperature 
dependence of the switching distribution width for the 
shunted nanowires. We find that it shows a trend similar 
to that previously observed for the unshunted wire[Io|. 
As for the retrapping current, which is observed to be 
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stochastic only in the shunted nanowire, the standard 
deviation is never observed to increase with decreasing 
temperature. However, at low temperatures, MQT is ex- 
pected to cause the standard deviation of the retrapping 
current to be constant with temperature [59]. In our ex- 
periments, we have seen evidence of this behavior and 
will investigate this more fully in the future. 



IV. THEORETICAL ANALYSIS AND 
DISCUSSION 

In this section we will develop a theoretical interpre- 
tation and understanding of the experimental data pre- 
sented in the previous section. In doing so, we will start 
by giving a brief account of the unshunted nanowires 
which have been previously investigated in detail. Then 
we will argue how shunting the nanowires brings about 
qualitative changes in the dynamics and develop a physi- 
cal picture by gathering experimental signatures and the- 
oretical arguments. Next we will motivate the theoretical 
model for explicit calculations and numerical simulations 
and use it to generate the V-I characteristics and the 
distributions that will be compared with experimental 
measurements. 



A. Unshunted nanowires 

Recently, properties such as the V-I characteris- 
tics and the switching distributions of the unshunted 
nanowire have been studied in detail to understand 
the behavior of quasi-one-dimensional superconductors 
at low temperatures [l(| El, US [Hj]. In quasi-one- 
dimensional superconductors the zero resistance super- 
conducting state is destabilized by thermal and quan- 
tum phase-slip fluctuations. These phase-slip fluctua- 
tions induce resistance which causes Joule heating in 
the nanowire. If this heat generated by phase-slip fluc- 
tuations in the bulk of the wire is not overcome suffi- 
ciently rapidly, it can reduce the depairing current to 
below the applied current, thus causing transition to 
the highly resistive state. It has been found in experi- 
ments with unshunted nanowire that while the distribu- 
tion of re-trapping currents is very narrow and almost 
temperature-independent, the distribution of switching 
currents is relatively broad and the mean as well as the 
width of the distribution change with temperature of the 
leads. The distribution in switching currents reflects that 
the collective dynamics of the superconducting conden- 
sate evolves stochastically in time and undergoes phase 
slip events at random instants. A stochastic model for 
the time-evolution of the temperature in a nanowire has 
been developed to understand above experimental results 
(27l l28j . The model predicts that although, in general, 
switching from superconducting to resistive normal state 
occurs due to several phase-slip events, it can be even 
induced by a single phase slip at a particular tempera- 



ture and current-range. The model also indicates non- 
monotonic temperature dependence of the width of the 
distribution of switching currents. Thus, these experi- 
ments with switching events as well as those with mi- 
crowave radiation in unshunted wires suggest that the 
resistive state is the normal state of the wire maintained 
by Joule heating, i.e. the JHNS [M EH, H HI • For un- 
shunted wires, the retrapping process is non-stochastic 
since the retrapping occurs from the thermalized Joule 
heating state. 



B. Qualitative picture of shunt-induced crossover 
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FIG. 7: Log-linear plot of average power dissipated in the 
nanowire immediately after a switching event (crosses) and 
just before retrapping (circles), plotted versus the shunting 
conductance (= 1/Rs)- Above a critical power, P c *, the resis- 
tive state of the nanowire is identified as the JHNS while be- 
low P* , the resistive state originates from a coherent phase dy- 
namics. The power is calculated by taking the voltage across 
the wire times the current through the wire (P* — IwV) at 
the retrapping and switching current in the resistive branch 
of the V-I curve from Fig[2ja). 

How does the picture for the unshunted case discussed 
above evolve as we shunt the wire with an external resis- 
tance? With the inclusion of a shunt resistance the ap- 
plied bias current is divided into two parts, and the part 
going through the nanowire decreases as the shunt resis- 
tance becomes smaller. A lower (higher) current through 
the wire causes a decrease (increase) in Joule heating in 
the wire and hence a decrease (increase) in local temper- 
ature. For the unshunted case that corresponds to an 
infinite shunt resistance, the heating is maximum and as 
discussed in the previous subsection one gets a JHNS in 
the unshunted wire. 

We will argue that upon decreasing the value of the 
shunt resistance the nature of the resistive state of the 
nanowire changes from a JHNS to a phase slip center 
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(PSC) state. PSC is a process of periodic-in-time phase 
rotation occurring in a certain region of the wire and is 
driven by the bias current. An ideal PSC acts qualita- 
tively like a weak-link JJ in series with the rest of the 
wire, and the differential resistance associated with it 
is determined by the quasiparticle diffusion length [olj . 
This resistance is much smaller than the normal resis- 
tance of the wire. This is what we indeed obtain when 
the nanowires are shunted with small resistance. The 
most explicit proof for a PSC would of course be the ob- 
servation of Shapiro steps, because they prove that there 
is a periodic phase rotation in the system. However we 
have compelling arguments and consistency in our expla- 
nation that points to the existence of PSC. 

The deterministic retrapping current observed for the 
unshunted nanowires becomes stochastic when an exter- 
nal resistive shunt is added (at least within the exper- 
imental range of shunt resistances). The deterministic 
retrapping current reflects that the resistive state of the 
unshunted wires is a thermalized JHNS. There is over- 
heating and the wire is normal. As I is reduced the tem- 
perature goes down. Relative fluctuations of the temper- 
ature are small since it is determined by a macroscopic 
number of normal electrons. In this case / r is fixed by 
the current at which the heating is not enough to keep 
the system above T c . If the system must re-trap at 1° 
then as can be seen in FigfJJ the distribution is sym- 
metric and Gaussian, centered around I®, and as wide 
as the noise in the bias current circuit and in the mea- 
surement circuit can smear it. On the other hand, the 
stochastic retrapping current indicates that the finitc- 
voltage/resistive state of the shunted wires is governed 
by a coherent dynamics of the phase of superconduct- 
ing order parameter. The dynamic state, called PSC 
moves by inertia, which is the voltage on the electrodes. 
But a strong thermal (or quantum) fluctuation can re- 
trap the system from the dynamic to the static state. 
And such change will be permanent. So the retrapping 
can happen at I > I®. The distribution is asymmetric 
(right-skewed or right-tailed) since the system can never 
switch at I < I®. Here, the fluctuation- free retrapping 
current Jj? is the current value at which the friction must 
stop the dynamic state in all cases (this is the property 
of the model considered). Similarly, the fluctuations in 
the low-resistance 'superconducting' state as discussed in 
the previous subsection for the unshunted case can allow 
for premature switching at < I < I® w = I c q but never 
at I > J,? w since the system cannot move through 7f w 
without a switch (the bias current I is assumed to grow 
linearly in time). One again has an asymmetric distri- 
bution, this time left-skewed or left-tailed. Overall, we 
can use the shape of the distribution to gain insights into 
nature of the state from which retrapping or switching 
happens. 

In Fig 13 a phase diagram for sample SI is presented 
which demonstrates the conditions necessary for the re- 
sistive state to be either the JHNS or a coherent phase 
dynamics state such as a PSC. The power P* = I\vV 



at switching and retrapping is calculated by taking the 
product of the current through the wire and voltage 
across the wire at which the system exhibits switching 
and retrapping, respectively. Here, 7w — I ~ ^ > where 
I is the total current, which obeys Kirchhoff's Law for 
current conservation (I = I s hunt + Iw> where /shunt is the 
current through the shunt). The critical power, P* , is 
defined as the minimum power the wire can sustain and 
still remain in the JHNS and is calculated from the power 
that the unshunted nanowire exhibits at retrapping. For 
sample SI, P* is calculated to be 0.533 nW from the un- 
shunted curve in FigJ5Ja). At switching, the unshunted 
wire experiences 31 nW of heating, which puts it in the 
JHNS, where it remains until the current is reduced be- 
low the retrapping current. 

With an included shunt, the Joule heating power in the 
wire at switching is reduced. For example, when shunted 
with 75 f2, the Joule heating power at switching is 0.359 
nW (compared to 31 nW for the unshunted wire), which 
is lower than P*. Thus, the wire switches to the PSC, 
which is a superconducting dynamic state (and not the 
normal state), and as the current is reduced, it remains 
in it until retrapping occurs. Because retrapping occurs 
from the phase coherent state, stochastic retrapping is 
expected for sample SI when shunted with (from the ex- 
perimentally examined values) a 75 il resistor or less. 
Some heating is also to be expected since the power at 
retrapping is still comparable to P* . 

Guided by the qualitative picture developed above, in 
the next subsection we will discuss a model that we will 
use for simulating the dynamics of the nanowire. In Sub- 
sec. IIVDI we will address the kinks seen in FigJSJa) and 
argue that they further support the existence of a coher- 
ent phase dynamic state or a PSC for shunted nanowires. 

C. Theoretical model for shunted nanowires 




FIG. 8: Circuit of the resistively capacitively shunted 
nanowire with an external shunt resistance Rs and the 
Johnson-Nyquist current-noises I n , 

As argued in the previous subsection, shunting changes 
the high-resistance state from JHNS to a state with co- 
herent phase dynamics or PSC. When the value of the 
shunt resistance is small, we expect that the dynamics 
of the nanowire can be modelled by an effectively zero- 
dimensional circuit model as is done in the case of a JJ 
given that the PSC behaves like a JJ in series with the 
wire. For unshunted nanowires (i.e. infinite shunt value) 
heating is important and as discussed above the dynamics 
of the nanowire is dictated by a model stemming from a 
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heat diffusion equation. For very large shunt values some 
elements of this model might still be important. However, 
the point of view we take is to simulate the dynamics of 
shunted nanowires using a circuit-element representation 
and see how well we can reproduce the experimentally 
observed behavior. 

As discussed in the introduction, the RCSJ model of 
Stewart and McCumber has been greatly successful in 
understanding the physics of JJs. So our goal is to adopt 
and extend this model to reflect the experimental set-up 
and measurements for the nanowires considered in this 
paper. The model of Stewart and McCumber was origi- 
nally introduced for superconducting JJs to study dc V- 
I curves displaying hysteresis for light damping [48l |49j . 
This model considered only the time varying phase dif- 
ference 4>(t) of the superconducting wave- functions in the 
weakly coupled superconductors and neglected any spa- 
tial variations of the superconducting wave-functions and 
is essentially zero-dimensional. We add an external shunt 
resistance in parallel with the superconducting junction 
in the RCSJ model to be able to study the behavior for 
different values of external shunt resistance as has been 
studied experimentally. We then simulate the extended 
RCSJ model with the Johnson-Nyquist Gaussian white 
thermal noise coming from the resistive parts of the cir- 
cuit. We have not included external noise in our sim- 
ulation. Before presenting the details of the model, we 
pause to provide some examples of analogies of the ob- 
served nanowire behavior with the established behavior 
in J J to further motivate the use of a RSC J kind of model 
for the case of a shunted wire. 

The hysteresis in underdamped JJs is due to the bista- 
bility of the phase point in the tilted wash-board poten- 
tial which depends nonlinearly on the bias current [3lj |. 
Damping plays an important role in dictating the dynam- 
ics of the J J. The experimentally observed saturation of 
isw at low shunt values as presented in Section. Mil can 
be interpreted to be an effect of high damping (damping 
~ Q^ 1 ~ Rg 1 ) on the premature switching process, and 
it indicates that the depairing current is nearly reached 
for these low values of the shunt. Another experimen- 
tal observation we presented earlier is that below some 
critical shunt value, the retrapping and switching current 
become equal and the hysteresis vanishes. For instance, 
for sample SI at 1.8 K at a shunt value of 10 f2, the V-I 
curve becomes non-hysteretic as in Fig. 2(a), and there 
is no abrupt switch into the resistive state. This can be 
interpreted in analogy to J.J as follows: At some critical 
value of the shunt, the increased damping changes the 
system from an underdamped junction (with hysteresis) 
to an overdamped junction (without hysteresis). In JJs, 
this transition occurs when Q ss 0.84 [38|. the third ex- 
ample of analogy with JJs is that the mean value of the 
retrapping current changes considerably upon changing 
the value of the shunt resistance. Indeed it is well-known 
for JJs that retrapping is very sensitive to the value of 
damping and the fluctuation free retrapping current is 
inversely proportional to the resistance associated with 



the J J. 

We give a brief summary of the physics of the RCSJ 
model in Appendix lAl and focus on discussing the details 
of our extended RCSJ model below. The displacement 
current and "normal" losses (e.g. quasiparticle tunnel 
currents) in the nanowire are included in the model by 
the shunting capacitance C and resistance i?w 5 respec- 
tively. We also include a Johnson-Nyquist type Gaussian 
white noise current source I n associated with the resis- 
tance i?w along with the drive current source I [H, [63[ . 
Next we extend the RCSJ model with an external normal 
resistance Rs and corresponding Johnson-Nyquist Gaus- 
sian white current noise I n ± for the present experimental 
set-up of a shunted nanowire (see FigJS]). Then, the re- 
duced equation of motion for the phase difference is given 
by (check Appendix [X] for a derivation) 



i + in + inl , 



(4.1) 



where Qo = (2eI c0 R 2 y C/h) 1 / 2 is the quality factor, 
i = I/Iqo is the normalized dc bias current, i a = 
In/Ico, W — Ini/IcO are the normalized noise cur- 
rents, t' = (2e/ c oi?w /K)t where t is the physical time 
associated with the circuit in Fig{8j Here J c o is the 
fluctuation-free critical current of the nanowire. The 
time-averaged steady-state voltage across the wire, V — 
IcoRw(d*4>/dt')t>, and the noise autocorrelations are 



(i n (t[)) = 0, (i nl (t' 1 ))=0 



hi 



cl) 



hi 



cO 



Rs 



(4.2) 
(4.3) 



where ( ) denotes averaging over the noise realizations 
(noise ensemble). The temperature of the nanowire and 
the shunt resistance are T\y and Ts, respectively. It is 
possible that the temperature of the wire is different from 
that of the shunt resistance, so we keep here two different 
noises coming from two different resistances. The rela- 
tions in Eqs. (|4.2l4.3[) are known as fluctuation-dissipation 
relations. Now if we assume that there is not significant 
MQT at the temperature (1.8 K) where the distributions 
for the switching current 7 SW and the retrapping current 
I r are measured, then the distributions are due to the 
thermal fluctuations arising from the Johnson-Nyquist 
current noises associated with the resistive parts of the 
circuit. 

One can write the Eq. (|4.1[) in a little different form as 



-rrj2 + ^7 + sin^ = i + i n + i nl , 



dt 



(4.4) 



where Q = (2eI cQ R^C/h) 1 / 2 , t' = (2eI c0 R T /h)t and 
1 /Rt = 1/ Rw + 1 /Rs ar e respectively the quality factor 
and the resistance of the full circuit. But the above form 
of Eq. (|4.ip is more useful in simulations. 
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D. Shunt dependence of V-I characteristics 




l(uA) 

FIG. 9: Voltage-current characteristics of the nanowire with 
external shunt resistance and current noise. The quality fac- 
tor Qo = 7 and the temperatures are related by Tw = Ts = 
T ~ 1.8 K and constant I c o = 5.55 fiA. The resistance of the 
unshunted nano-wire R n = 1385 fl, and the resistance of the 
phase-slip center (or shunted nano-wire) -Rpsc ~ 0-lR n - 

Now we simulate Eq. (|4.1[) for T\y = Ts = T e i ectro de 
to calculate the voltage-current characteristics of the 
shunted nanowire in the presence of the current noises 
from the normal resistances. In the experiment, one 
changes the bias current with a finite current sweep rate, 
and measures the corresponding voltages. In the simu- 
lation, instead we first fix a bias current and then in- 
tegrate the above equations of motion (with suitable 
initial conditions and current noises) for a sufficiently 
long time (this time is the relaxation time or the tran- 
sient time), and next calculate the time-averaged volt- 
age by averaging d<p/dt' over some time interval. For 
forward current sweep we choose the initial conditions, 
(j)(f = 0) = 0,^(f = 0) = 0; and for the backward 

current sweep we use <j>{if = 0) ^ 0, J$(£' = 0) ^ 0. 
The initial values of 4>(t') and J^(i') for the backward 
current sweep can be any finite non-zero values in the 
resistive state of the system as after long time of tran- 
sient dynamics the exact initial values of 4>(t') and tt(t') 
are irrelevant. We generate the Gaussian white noises i n 
and i n i at each time step of the simulation satisfying the 
noise properties in Eqs. (|4.2[ 14. 3|) following the method 
described in [64| . 

The voltage-current characteristics are plotted in FigJ9] 
for the shunted wires with different values of the normal 
resistance of the wire and the shunt resistance. As can 
be seen from the Eqs. (|4.1IA5[) , we really do not need ex- 
plicit values of the resistances in simulation, but we only 
need the ratios of the two resistances and a quality factor. 
We have checked that for the quality factor Q = 7 and 
at a specific temperature 1.8 K, the transition from hys- 
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FIG. 10: Comparison of the voltage-current characteristics of 
the shunted nanowire with Rs = 10 Q from the experiment 
and simulation. We find the best agreement from the sim- 
ulation of the shunted nanowire for a quality factor Qo = 7 
and temperature Tw = Xs = T ~ 0.13 K, I c o = 5.15 /J.A, 
and a resistance of phase-slip center (or shunted nanowire) 

r psc = 85 n. 

teretic to non-hysteretic behavior occurs near the ratio 
i?w : i?s = 4 : 1. As discussed in the previous subsec- 
tion, our present theory with a coherent phase relation- 
ship suits us best to analyze the nanowires shunted with 
low shunt resistance, i.e., in the PSC regime. A compar- 
ison between FigJ5Ja) and Fig |H] shows a very good qual- 
itative agreement between the experimental V-I curves 
and those from simulations. We also have good quan- 
titative agreement between the experiment and simula- 
tion for the V-I curves of the shunted nanowires (PSC 
regime) with small resistances. This is shown in Fig llOl 
for the shunted nanowire case where Rs = 10 f2. We 
find the best agreement between the experimental and 
simulated V-I curves for the low shunt resistances when 
the effective temperature (related to the thermal noise) of 
the simulated nanowires is much lower than the tempera- 
ture of the superconducting electrodes in the experiment. 
Therefore the fluctuation in 7 SW due to thermal noise is 
expected to be small as we find in experiment. The rea- 
son for such effective thermal noise reduction is not com- 
pletely clear. It might be due to the inductance of the 
shunt resistance. The inductance can cut-off some of the 
higher-frequency thermal noise, thus reducing the stan- 
dard deviation of the noise current. It also confirms that 
Joule heating in the shunted nanowires for lower shunt 
resistance is greatly reduced compared to the unshunted 
case. The resistance of the shunted nanowires i?w which 
enters in the simulation through Eqs l4.ll 14.31 is used as a 
fitting parameter here. It is necessary to choose it to be 
much smaller than the normal resistance of the nanowire 
R n to have the best fitted of V-I curves. This is a strong 
indication that the shunt resistance drives the nanowire 
to a phase-coherent PSC state, in which the time-average 
supercurrent is not much smaller than the total bias cur- 
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rent. Thus we introduce a new notation for the wire 
resistance, namely -Rpsc- This quantity represents the 
value of the wire resistance that we have to put in our 
model in order to produce the best fits to the experimen- 
tal V-I curves. We find that i?psc << Rn- For example, 
we find i?psc = 85 from the simulation for Rs = 10 fl 
case. Here we remind that the normal resistance of the 
nanowire is R n = 1385 Q. 

In the experimental V-I curve of the shunted nanowire 
with shunt resistor 25 f2 there are kinks which we do not 
find in simulations. The kinks can be attributed to the 
effects of a shunt inductance in series with the shunt re- 
sistor. These kinks are not associated with resonance in 
the system because such a resonance would not depend 
on temperature as these do. Such inductive effects orig- 
inate from the fact that the resistor used for shunting 
has dimensions of a few centimeters and so has a large 
inductance (~ 20 nH). Inductance connected in series 
with a shunt resistor is known to cause similar kinks in 
the V-I curves of shunted JJs due to a complicated dy- 
namic of the phase difference on the junction [65}. Thus, 
the observation of such kinks confirms that the resistive 
state in our shunted wires is due to a phase-coherent PSC 
[3lj and not due to Joule heating. Thus we find another 
indication that by resistively shunting the nanowire it is 
possible to change the nature of its resistive state from a 
phase-incoherent JHNS to a phase-coherent PSC state. 

E. Shunt dependence of switching and retrapping 
distributions 
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FIG. 11: Simulated distribution P(7 SW ) of the switching cur- 
rent 7sw for different Rs/Rw values at temperature Tw = 
Ts = 1.8 K, constant To = 5.55 /xA. The quality factor 
Qo = 7. Inset: Shows the standard deviation <t sw of switching 
current distribution as a function of Rs/Rw- 

We first calculate the distributions P(/ sw ) 01 the 
switching current (I aw ) in the shunted nanowire for the 
same temperature Tw = Tg — T e i ec trode- We again simu- 
late Eq. (|4.ip as before by changing the bias current with 



zero as initial values for <fi and , but now we repeat the 
full procedure for many realizations of the thermal noise. 
We count a switch from the metastable to the running 
state when the nanowire spends more than half the time 
in the running state over some sufficiently long time pe- 
riod (t s ) of simulation. This gives us a distribution for I sw 
which is plotted in Fig QT] for the shunted wire. We find a 
good qualitative agreement between the simulation and 
experiment for the switching distributions of the shunted 
nanowires as shown in FigJ3] and Fig llll In the inset of 
Fig llll we show the standard deviation <r sw of the simu- 
lated P(I SW ) with Rs/Rw and it matches with the trend 
of the standard deviation of the measured P(J SW ) from 
the experiment. Here we mention that cr sw shows a non- 
monotonic behavior with Rs /Rw in the RCS J model for 
higher values of Rs, for example, Rs > Rw at a constant 
temperature. This non-monotonic behavior of er sw in the 
RCS J model with Rs/Rw at a constant temperature is 
similar to non-monotonic behavior of <r sw with tempera- 
ture for a constant Rs/Rw- It will be discussed in the 
last subsection. However, experimental results show the 
value of <t sw is greater for the unshunted case (Rs = oo) 
than the shunted case. This also indicates that the un- 
shunted nanowire is not in a coherent PSC state but is 
dominated by Joule heating which increases the effective 
temperature of the wire. 
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FIG. 12: Simulated distribution P(I r ) of the retrapping cur- 
rent I r for different Rs/Rw values at temperatures Tw = 
Ts = 1.8 K, constant I c o — 5.55 fj,A. The quality factor 
Qo — 7. Inset: Shows the standard deviation cr r of the re- 
trapping current distribution as a function of Rs/Rw- 

We next simulate the extended RCSJ model to under- 
stand the measured distributions P(I r ) of the retrapping 
current for different shunt resistances. Here we choose a 
fixed temperature. The simulation method is similar to 
finding the switching distributions, but now we start from 
a PSC state with nonzero initial conditions for <p(0) and 
3^(0). We reduce the bias current and count a retrap- 
ping event from the running to metastable state when 
the nanowire spends less than half of the time in the 
running state over time period of r s . The simulated re- 
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trapping distributions are plotted in Fig[12] for different 
Rs/Rw values at a constant temperature. We find from 
Fig|4]that the standard deviation of the retrapping cur- 
rent distributions falls slightly (within the experimental 
noise limit) with decreasing shunt resistance. The stan- 
dard deviation of the simulated retrapping current distri- 
bution also falls slightly below Rs/Rw ~ 0.4 for Q = 7 
which is similar to the experiment. But we also find 
from the simulation of the RCSJ model that the stan- 
dard deviation of P(I T ) decays slightly with increasing 
shunt resistance above Rs/Rw ~ 0.4 for Q — 7, but 
it never goes to zero at higher shunt resistances for the 
RCSJ model with coherent dynamics. The standard de- 
viation of the measured retrapping current distribution 
for the unshunted wire is almost zero within the noise 
limit; this is consistent with the existence of JHNS in the 
unshunted nanowire and PSCs in shunted wires. The 
width of the simulated switching and retrapping distri- 
butions are much greater than the experimental results 
at the same temperature T = 1.8 K. This might be again 
due to the inductance of the shunt resistor which can 
effectively reduce the thermal noise in the system. 



F. Temperature dependence: shunted vs. 
unshunted nanowires 



tuations has also been obtained previously for various 
JJs [36l - l38j and our theoretical analysis not only high- 
lights its ubiquity but also provides a way of obtaining 
it in terms of a RCSJ kind of modelling. We also find a 
non-monotonic temperature dependence of the standard 
deviation of retrapping current in our simulation. In our 
numerical study here we only consider effect of thermal 
fluctuations in phase-slips, thus one needs to go beyond 
the present study to include effect of macroscopic quan- 
tum tunneling on the temperature dependence of ct sw in 
the fully quantum regime. 

In unshunted nanowires too it has been shown that the 
standard deviation of the switching current distribution 
is non-monotonic as a function of temperature [27],[28j|. In 
the thermal regime at higher temperatures T > T* multi- 
ple phase slips are required before the wire switches to the 
normal state and the standard deviation increases as the 
temperature is decreased [l(| ■ At slightly lower tempera- 
tures T < T* a single thermal phase slip causes the wire 
to switch to the normal state and the standard deviation 
decreases with a decrease in temperature. At low tem- 
peratures when QPS are present, depending on how T* 
compares with the temperature of crossover from TAPS 
to the QPS, one can get different behaviors. With an ap- 
plied external shunt, the increased dissipation is expected 
to decrease the temperature at which the crossover from 
thermal activation to MQT takes place. 
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FIG. 13: The standard deviation <j sw vs. T for the simu- 
lated switching events in the RCSJ model with Qo — 7 and 
Rs/Rw = 0.3. 

Finally we simulate the extended RCSJ model with an 
external low shunt resistance to find the temperature de- 
pendence of the standard deviation er sw of switching cur- 
rent. We use the same scheme as previous sub-sections 
to determine the switching current in the simulation at 
different temperatures. We plot ct sw vs. T in FigfTS] 
for Rs/Rw — 0.3 and Qo = 7. We find that temper- 
ature dependence of simulated a sw due to the thermal 
fluctuations is non-monotonic just as in the experiment 
(see FigjS]) for the shunted nanowires. A non-monotonic 
temperature dependence of a sw due to the thermal fluc- 



V. CONCLUDING REMARKS 

We have undertaken a detailed study of the effect of ex- 
ternal resistance shunts on the behavior of superconduct- 
ing nanowires. Shunting has a strong effect on the be- 
havior of the nanowire. We find that the statistics of the 
switching and retrapping currents significantly depends 
on the value of the shunt resistance. The temperature 
dependence of the mean switching current in strongly 
shunted nanowires is consistent with the Bardeen predic- 
tion for the temperature dependence of the critical cur- 
rent [58| ; this indicates that the switching current can 
be controllably driven very near to the depairing cur- 
rent through external resistive shunting. The retrapping 
current on the other hand increases and becomes more 
stochastic, at least for moderate shunting. We demon- 
strate that the shunting, even with a large resistance 
value, can be used to control the phase slip events in 
the wire. We suggest a model based on the Stewart- 
McCumber RCSJ model, which is generalized to include 
two resistive elements, corresponding to the effective re- 
sistance of the wire (with a phase slip center), and the 
resistance of the shunt. The model provides a semi- 
quantitative description to the data. Moreover, it pro- 
vides insights into developing a circuit-element represen- 
tation of a superconducting nanowire. 

Our work opens up many interesting avenues towards 
developing a fundamental understanding and control of 
coherence and dissipation in nanowires as well as its rela- 
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tion to possible quantum phase transitions. It will be im- 
portant to develop a model that incorporates heating as 
well as coherent dynamics and studies the entire crossover 
in going from unshunted nanowires (i.e. infinite shunt re- 
sistance in parallel) where heating is most important to 
the case of very low shunt resistance values where heating 
is least important. It would be extremely interesting and 
relevant also to explore the low temperature quantum 
regime in more detail and to study the implications of 
our work for quantum computing and other technologi- 
cal applications of nanowires. As an example, nanowires 
have been used as photon counters [66j, which are im- 
portant in radioastronomy. The switching events studied 
in this paper, represent so-called "dark counts", in the 
terminology of the photon detection community. Under- 
standing the physics of dark counts is important for the 
purpose of improving superconducting photon detectors. 
The fact that the standard deviation of switching current 
becomes smaller with the inclusion of a shunt resistor has 
relevance to photon detectors, since dark counts can be 
reduced by shunting. 
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Appendix A: Resistively and capacitively shunted 
Josephson junction (RCSJ) model 




FIG. 14: Stewart-McCumber resistively and capacitively 
shunted Josephson junction (RCSJ) model with the Johnson- 
Nyquist thermal current noise. 

Here we briefly digress the main features of the RCSJ 
model (see Fig fl4]) introduced by Stewart and McCum- 



ber [48|, |49(. We find the equation of motion for the 
time varying phase difference <fi(t) of the superconduct- 
ing wave-functions of the circuit in Fig ll4l b y apply ing 
the well-known Josephson dc and ac relations [67H69( for 
the current-phase (I — </>) and the voltage-phase (V — </>), 



I 3 = IcO sin 4> , 
d6 _ 2eV 
~dt ~ ~h~ ' 



(Al) 
(A2) 



where I c o is the fluctuation-free intrinsic critical current 
of the junction. The equation of motion for <p{i) of the 
circuit in Fig ll4l is given by, 



Ch d 2 (j) h d(j) 
+ 2eRi~dl 



I cQ sin < 



/ + /n 



(A3) 



along with the Gaussian white noise properties for the 
Johnson-Nyquist thermal current noise J n , 



&(*)>= 

(/n(tl)Jn(*a)> = ^jpS(t! 

R.l 



(A4) 



where Xj is the temperature of the JJ with capacitance 
C and resistance Rj. The resistance Rj measures dis- 
sipation in the JJ in the finite voltage regime, without 
affecting the lossless dc zero voltage regime, and C in- 
dicates the geometric shunting capacitance between the 
two superconducting electrode [31|. The Eq. (|A3p of mo- 
tion for the junction phase can be rewritten in terms of 
dimensionless parameters as 



2 d 2 cj) d(j> 



'dt>' 



dt> 



l + in, 



(A5) 



where Qq = (2eI c oR 2 C/h) 1 ^ 2 is the quality factor of the 
linearized equation of motion, i = I/I c o, in — In/IcO, 
and t' = (2eI c oRj/h)t is the normalized time. The 
term t£ is damping as it breaks the time-reversibility 
of the equation and introduces dissipation. The strength 
of damping is proportional to l/i?j and is inversely re- 
lated to the quality factor. In this notation, the time- 
averaged steady-state voltage across the junction, V — 
IcoRj(d<fi/dt')t>- The noise correlation in the scaled time, 
(i n (t'iK(i 2 )) = (4efc B T w /nJ c o)5(ii - t' 2 ). The usual 
McCumber parameter f3 c = Q\ and Stewart parameter 
ijJqt = Qq with ujq = y/2eI c o/HC and r = i?jC For the 
circuit in Figj8] we replace i?j in Eq. (|A3|) by Rt where 
l/R T = l/R w + l/R s . Then we derive either Eq. (|4~Tj) or 
Eq. (|4.4p following the similar steps to get Eq. (|A5|) from 
Eq.dMl. 

In the absence of thermal current noise 7 n at zero 
temperature (also neglecting quantum fluctuations), the 
zero-voltage state or state is stable at all bias levels 
less than the ideal critical current (\i\ < 1), and the volt- 
age state or 1 state is stable at all bias levels greater 
than a minimum value designated by a fluctuation free 
re-trapping current z r0 . The value of i r o (= Iro/I c o) is de- 
termined entirely by the quality factor Qq and decreases 
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with increasing Qq as a smaller tilt is sufficient to support 
the running (finite voltage) state when damping is less. 
For Qq < 0.8382, the damping is sufficient that a running 
state is not possible unless the potential decreases mono- 
tonically and in this case i r o = 1. For Qo > 0.8382, a run- 
ning state is possible even when the potential has local 
minima. [70l ] In this case « r o < 1 and the V-I curve is hys- 
teretic. In the limit of large Qo, *r0 = 4/7rQ (Qo > 3). 
48, 70] 



The phase dynamics described in Eq. (|A5p can be vi- 
sualized as the damped motion of a Brownian particle in 
the tilted washboard potential U(4>) — —{i(j> + cos^). In 
the under-damped regime Qq > 1, the zero- voltage state 
and the resistive state correspond to the particle trapped 
by the energy barrier AU and running downward along 
the tilted potential, respectively. Escape from the poten- 
tial (0 state to 1 state) can occur even for i < 1 due to 
the thermal and the quantum fluctuations. 
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